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In this talk I review the low x QCD experimental results from HERA and their 



theoretical underpinnings. In particular, I discuss the structure function F2, large 
rapidity gap events, pomeron structure functions and J/ip production. 
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1 Introduction 

The electron (positron)-proton collider facility at DESY, nicknamed HERA (for Hadron 
Electron Ring Accelerator) has, since its inception in 1992, given many new and interesting 
results from the two experiments HI and ZEUS. These experiments have thrown up a 
wealth of new theoretical issues pertaining to low x QCD. These experiments cover both 
the deep-inelastic region (DIS) and the photoproduction region. The basic interaction 
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Figure 1: The basic DIS diagram with kinematic variables labelled 



picture is the usual DIS diagram shown in figure 1 for the process 



e(k)+p(p)^e(k')+X. 



The following kinematic variables completely describe the process: 



s = 




^ 4E e E p 



2E e E e ,(l + cos9) 



(1) 



(2) 
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Q 2 Q 2 

x = - — ~ — 4 
2p.q ys 

W = (q + p) 2 ~ -Q 2 + ys (5) 

Here s is the CM energy squared, E e (E' e ) is the energy of the initial (scattered) electron, 
9 the angle of the scattered electron in the lab frame with respect to the proton direction, 
Q 2 is the negative of the momentum transfer squared to the proton, y is the fraction of the 
electron energy carried by the virtual photon in the proton rest frame, the Bjorken variable 
x is the fraction of the proton energy carried by the struck quark and W is the hadronic 
invariant mass. The masses of all particles are neglected. 

At HERA, a 27 GeV electron beam collides head-on with an 820 GeV proton beam 
giving 

s ~ AE e E p ~ 10 5 GeV 2 

which is much larger than hitherto obtained at fixed target experiments. (Since 1994, the 
electron beam has been replaced by a positron beam at the same energy). As a result the 
two experiments HI and ZEUS can measure the structure function F 2 in a completely new 
(x, Q 2 ) domain (see Figure 2). 

The Born cross section for single photon exchange in DIS is given by 

d 2 a 2ira 2 



dxdQ 2 Q 2 x 



V 2 

2(1 -y) + 



F 2 (x,Q 2 ) (6) 



1 + R_ 

where R is the photoproduction cross section ratio for longitudinal and transverse polarised 
photon, R = olIot- R has not been measured at HERA; therefore a QCD prescription 
based on parametrisation of parton densities is used. Q 2 is directly measured from the scat- 
tered electron but x is calculated from Q 2 and y and therefore depends on the experimental 
resolution of y. 



2 Measurement of F$ (x, Q 2 ) 

The data for Ff from HI and ZEUS |l], @ is shown in figures 3, 4, and 5, along with 
some of the older fixed target data. It is clear that at fixed Q 2 , F 2 rises with decreasing x 
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Figure 2: Phase space coverage of the F 2 measurements (from |]T3| ) 

down to the smallest Q 2 , the steepness of the rise decreasing with decreasing Q 2 . Similarly 
at fixed x < 0.1 F 2 rises with Q 2 , the rise becoming steeper as x decreases. (We shall 
henceforth drop the superscript p on F 2 since all our discussion henceforth will refer to 
the structure function of the proton). In fact, as seen in figure 4, upto Q 2 ~ 0.85Gey 2 
the data favour a parametrisation based on a soft pomeron such as that suggested by 
Donnachie and Landshoff J|. For larger values of Q 2 (Q 2 > l.OGeV 2 ) the usual QCD 
based parametrisations of GRV or MRS M, Q give a good description of the data. 

The behaviour of F 2 as a function of Q 2 for selected values of x is presented in figure 5. 
The growth of the structure function with Q 2 for x < 0.01 over three orders of magnitude 
in Q 2 is clearly seen. With the new data from HERA the gap between fixed target and 
the HERA experiments has been filled and the remarkable agreement with a NLO QCD fit 
performed by HI over three decades in Q 2 and over four decades in x is clearly seen. 

The HI experiment has parametrised the rise of F 2 with x by F ~ x~ x at fixed Q 2 
values. The result for A as a function of Q 2 is shown in Figure 6. It is clear from the 
figure that A increases with increasing Q 2 reaching a value of around 0.4 — 0.5 around 
Q 2 ~ 10 2 - 10 3 GeV 2 . 
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Figure 3: F 2 (x, Q 2 ) as a function of x. The curves are the NLO QCD fit for Q 2 > 5GeV 2 . 
Curves below 5 GeV 2 are obtained by backward evolution, (from 



3 The rise of F<i and perturbative QCD 

The structure function F 2 {x, Q 2 ) is given in terms of the parton distributions 



F 2 (x,Q 2 



x 



= J2 e t c i ® (ft + ft) + C g ®g. 



(7) 



The evolution of the structure function F 2 as a function of Q 2 is governed, in perturbative 
QCD, by the Dokshitser-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation || |7j, 
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Figure 4: F 2 (x,Q 2 ) as a function of x for low Q 2 compared with model predictions. (from 
0) 

Here q and g are the quark and gluon distributions, <S> denotes convolution with respect to 
x, i.e.[/ ® g](x) = $1 -£f(y)g(y), n f is the number of flavours of quarks, e, is the electric 
charge of the quark gj. The coefficient functions C's at leading order are 

Q(x,Q 2 ) =5(l-x); C g (x,Q 2 ) = 0. (9) 

At higher orders, they depend on the specific factorization scheme. In the common parton 
scheme the above remains true to all orders. 
The splitting functions have the expansion 

P l3 {x,Q 2 ) = %P§\x) + (^) 2 4 2) W + • • • • (10) 
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Figure 5: Compilation of Fi measurements as a function Q 2 for selected values of x. The 
numbers in parenthesis are the scaling factors by which the value of F2 has been multiplied 



in the plot. The curves are a NLO QCD fit by HI. (from |13 



The first two terms above define the NLO DGLAP evolution. 

The rise of F2 as a function of x is explained by the usual DGLAP evolution and also 
by more non-conventional dynamics like the BFKL evolution ||. Let us briefly review the 
two different approaches. 

To understand the leading log summation technique in the DGLAP approach, we use 
the Dokshitzer method 0. At low values of x, the gluon is the dominant parton in the 
proton. He showed that the GLAP summation of the (a s logQ 2 ) n terms amounts to a sum 
of gluon ladder diagrams of the type shown in Figure 7 with n gluon rungs. The form of the 
gluon splitting function P gg ~ 6/x for small x gives the evolution of the gluon distribution 
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Figure 6: Variation of the exponent A from fits of the form F 2 ~ x A at fixed Q 2 values 
and x < 0.1. (from 



function as 




It is clear that the (logQ 2 ) n builds up from the nested k integrations. In fact this con- 
tribution comes from a region where the transverse momentum of the emitted gluons are 
strongly ordered 

Q 2 » k 2 nT » ...» k 2 2T » k\ T (12) 

Similarly, the (log(l/x)) n term comes from the region where the longitudinal momentum 
fractions are strongly ordered 

X «£„«...«& (13) 
In the region of small x the log(l/x) terms have also to be summed. With the assumption 
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Figure 7: A gluon ladder diagram that contributes to the DLL or BFKL ladder summations. 



of £i<?(£i,(5o) approaching a constant, say G the result we get is 



xg(x,Q 2 ) 



'q_ 2 ' 



G exp < 2 



3a, . /Q 2 \ . 
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(14) 



in the limit of large Q 2 and small a;. 

This is an example of an all order double leading log (DLL) summation of a s log Q 2 log(l/x) 
obtained by summing the strongly ordered gluon ladder diagrams as we have just explained. 
It tells us that as x decreases, xg(x,Q 2 ) increases faster than any power of log(l/x) but 
slower than any power of 1/x. 

If kiT ~ k i+ iT then we lose one power of logQ 2 . However, at small x, this can be 
compensated by a large log(l/x) factor. Relaxing the strong ordering constraint on the 
/c T 's and summing just the log(l/a;) in the small x region will give us the leading log 
summation in log(l/:r) instead of DLL and we get 



xg(x,Q 2 ) 



V 1 7i 1 n\ 



c log - 

X 
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-X 
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x~ x (15) 



where A = (3a s /7r)c. This difficult summation was done rigorously by Balitskij, Fadin, 
Kuraev and Lipatov (BFKL) ||. The constant c = 4 log 2 and so, for a s ~ 0.2, A ~ 0.5. 
We have highly simplified the discussion of the BFKL equation, and in practise one needs 
to work in terms of the unintegrated gluon distribution f(x, k%) and integrate f(x, k\}fk\ 
to get the final gluon distribution. 

We thus have two different predictions for the rise of F2, the DLL (from the DGLAP 
evolution equations) and the BFKL. In principle it should therefore be possible to distin- 
guish between these predictions by looking at the F2 data from HERA. Unfortunately this 
is not that simple. The steepness of the rise of F 2 with decreasing x can be controlled by 
varying Ql, or the starting distribution <?(£i, Ql). In addition, the set of equations based 
only on F2 measurements is underconstrained and one needs one other measurement like 
the longitudinal structure function Fl in order to constrain the system fully. 

The net consequence of the above is that the present data on F2 does not distinguish 
betweem the BFKL and DGLAP (or for that matter the CCFM which embodies both) 
predictions. 

We would like to mention here a suggestion by Mueller || on measuring an observable 
that is less inclusive than the F2 measurement (in which none of the properties of the 
hadronic final state is measured). He suggests studying DIS events at small x containing 
an identified high-energy jet emitted close to the jet of proton fragments (see Figure 8). 
The identified jet originates from the parton labelled a. If we study events with Xj large 
(> 0.05) and x very small (~ 10~ 4 ) the ratio x/xj will be sufficiently small to reveal the 
(x/xj)~ x behaviour of the BKFL ladder. Details may be found in the references mentioned 
above. This study is currently under way and preliminary results are reported in JIT]]. 

The behaviour of F 2 as a function of both Q 2 and x (for large Q 2 and small x) can 
be neatly combined through the double asymptotic scaling analysis of Ball and Forte |[L1|| . 
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Figure 8: A ladder diagram contributing to an identified forward jet as discussed in the 
text 



They use the two variables 



a = \/tog(xo/x) log(a 8 (Q )/a 8 (Q)), p = 



\og(x /x) 



\ \og(a s (Q )/a s (Q)) 

with a s (Q) the strong coupling constant evaluated upto two loops. 

The parameters xq and Qq have to be determined experimentally. Qq is found to be 
2.5GeV 2 and x = 0.1. (See, for example F 2 is rescaled by the leading asymptotic 
factors 

R F (a, p) = Z.lexp ( -2 7 a + u- + \ logoff) + log(^) ] /f F (17) 
V P 2 7 / 

and 

R' F (a,p) = R F exp(2^a) (18) 

where 

i F = 1 + ((6 + 6H(Q) - ei«.(Qo)) W(2tt7))- (19) 
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The above constants are £1 = (206n//27 + 6bi/b )/bo, £2 = 13, b = 11 — 2n//3, 



(ll + 2ri//27)/6o and 61 = 102 — 38n//3. Also, n/ is the number of flavours and 7 = yl2/&o- 

Thus the function \og(R' F F2) is predicted to rise linearly with a and RfF2 is expected 
to be independent of p and a. 

Figure 9 shows the data for these two predictions. The first graph is for Q 2 > 3.5GeV 2 . 
Approximate scaling is observed for Q > 5GeV 2 and p > 2. At high p and low Q 2 , the data 
tends to violate the predicted scaling behaviour. The second graph shows the appropriate 
data for p > 2 and Q 2 > 5GeV 2 as a function of a. The predicted linear rise is seen with 
slope of 2.50±0.02±0.06 in good agreement with the QCD expectation of 2.4 for 4 flavours. 

Thus in the low x , high Q 2 region beyond Q 2 > 5GeV 2 , scaling is seen in the two 
variables p and a. The prediction of double scaling is in very good agreement in this 
region with experiment. Beyond the double scaling region, the evolution continues to be 
in excellent agreement with the full NLO QCD evolution prediction. 

One of the by-products of the measurement of the structure function F2 and its log Q 2 
evolution is that the gluon distribution inside the proton has been measured with greater 



accuracy (cf. Figure 5 in fT3|j). 

To summarise this section, the structure function F 2 has been found to rise steeply 
with decreasing x down to Q 2 as low as 0.85GeV 2 . This rise can be parametrised by a 
form x~ x where A is found to be a function of Q 2 and reaches a value of around 0.4 — 
0.5 around lOOOGel^ 2 . This behaviour is described very well by conventional NLO QCD 
(through the DGLAP equation). However this does not rule out the existence of non- 
conventional dynamics like BFKL because the system is underconstrained. Thus some 
other measurement which is less inclusive than F2 like the forward jet mentioned above, 
or a measurement of the longitudinal structure function F^ is required. In any case, some 
new physics must take over at sufficiently low x otherwise the indefinite rise of F2 would 
violate the unitarity bound. A study of the longitudinal structure function in the double 
scaling limit may be found in ||12||. 
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Figure 9: The rescaled structure functions a) RpF 2 vs. p and b) log(i?^F 2 ) vs. a as 
explained in the text. Only data with Q 2 > 5GeV 2 and p > 2 are shown in b). (from [|TJ) 



4 Large Rapidity Gap Events 

In approximately 10% of the events used to study F 2 , a large rapidity gap in the hadronic 
final state was observed between the photon and the proton fragmentation regions. This is 
described in Figure 10. These events were fairly independent of W and Q 2 . This implies 
that there is no energy or colour flow between the struck parton and the proton remnant 
direction. A generic explanation of these events is through the exchange of a colour singlet 
parton which is typically a pomeron. Such events are typical of diffractive events through 
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which one introduces the idea of the pomeron. For our purposes, the pomeron is the generic 
name for a particle the exchange of which is responsible for these large rapidity gap events. 



Figure 10: Schematic representation of a DIS event with a large rapidity gap as discussed 
in the text 

One of the simplest ways to understand these processes is to consider j*p scattering in 
the rest frame of the proton. Here, the life time of a qq fluctuation of the virtual photon is 
given by r ~ l/(2m p x) where m p is the mass of the proton. This time scale corresponds 
to a distance of cr ~ 10 3 /m for x = 10~ 4 which is much larger that the typical size of a 
hadronic target. Thus the photon may fluctuate into a qq pair before arriving at the target, 
thereby making 7*p interaction similar to hadron hadron interaction. Diffractive scattering 
emerges as a consequence of the hadronic nature of the virtual photon. 

It is important here to present the appropriate kinematic variables for the process. 
We use the variables xp for the fraction of the proton momentum carried by the generic 
pomeron iP, and (5 the fraction of the pomeron momentum carried by the struck quark 




p 



P' 
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with f3 = x/xjp, x being the usual Bjorken variable. Neglecting the mass of the proton and 
taking t ~ 0, (see Figure 10) 

Q 2 + M x Q 2 
XF ~ W 2 + M 2 X ] P ~ Q 2 + M 2 X (20) 

Assuming the validity of Regge factorization for 7*p diffractive interactions we have 

^l> =FF/p{t , XFMY1 P), (21) 

where Fjp/ p (t,xp) is the flux of the pomeron in the proton. The DIS differential cross 
section for producing a diffractive hadronic state X of mass M x is 

^^^[l + d-^^.O.,,,), (22) 

denotes the contribution of diffractive events with a given xjp integrated over t to 
the proton structure function F2 for a given Q 2 and x. This part of the structure function 
is related to the pomeron structure function Ff by 

F 2 D{3 \p,Q 2 ,x P ) = -\-Ff(P, Q 2 ) , (23) 

Xjp 

where n = 1 — 2ajp(0) + 5 t independent of (3 and Q 2 . 

HI has carried out measurements of xjpF®^ in the range 2.5 < Q 2 < Q5GeV 2 , 0.01 < 
(3 < 0.9 and 10 -4 < xjp < 5.10 -2 and this is shown in Figure 11. The pomeron trajectory 
fit corresponds to an intercept 

ajp(0) = 1.18 ±0.02 ±0.07 (24) 

which is somewhat higher than the intercept of a soft pomeron (~ 0.08). The result from 
ZEUS is similar 

ajp(0) = 1.17 ± 0.04 ± 0.08 (25) 

The quantity F® = J F^ } dxp is directly proportional to F^ and is found to be almost 
independent of Q 2 . This can be explained only by a large gluon contribution (see for 



details). In fact more than 80% of the Pomeron momentum is carried away by hard gluons 
in the range of Q 2 being studied. 
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To summarise, the pomeron, which is a generic term to describe a colour singlet particle 
with vacuum quantum numbers, seems different in different processes. In particular, the 
pomeron in these diffractive events with large rapidity gaps seems to be harder than the 
usual soft pomeron of Donnachie and Landshoff ||. Using DGLAP evolution it is found 
that the effective pomeron is made of mostly hard gluons. It is clear however, that the 
pomeron is not a universal object and has different intercepts in different processes. 



5 Photon Structure Function 

Even though a photon, being a gauge particle has no intrinsic structure, its hadronic fluc- 
tuations (through its interaction with matter) give it a partonic structure. Hard scattering 
processess (jets and charm production) yield tests of perturbative QCD at NLO and give 
information about the quark and gluon structure of the photon. We will very briefly touch 
upon some of these matters. Large dijet photoproduction takes places at leading order 
through the 'direct' scattering process a(jq — > qg)®q p and a{jq — > qq)®g p . There are also 
so-called resolved processes (at higher order) a(gq — > gq)®g^®q p and &(qq — > qq)®q p ®q~ f . 
The momentum fraction of the photon constituent can be recontructed from the transverse 
energies and rapidities of the dijets, 

_ gri exp(-77i) + E T2 exp(-^ 2 ) 

Thus, 'direct' processes have x 7 ~ 1 and 'resolved' processes x 7 < 1. The angular dis- 
tribution of the dijets can also be measured to provide information on the nature of the 
hard scattering process. In fact the expected difference in distribution between direct and 
resolved process can be clearly seen in the data and is again in good agreement with NLO 
QCD predictions. There are also a host of other results both from HI and ZEUS. A good 
summary of the results and their theoretical implications can be found in the report of the 
photoproduction working group of the HERA workshop on proton, photon and pomeron 
structure function WA . 



15 



6 Exclusive vector meson production 



The detectors at HI and ZEUS have the ability to study the nature of the hadronic final 
state in detail. As a result there has been renewed interest in the study of J/0 production 
at HERA. The production cross section for the process is dominated by photoproduction 
i.e by the interaction of almost real photons Q 2 ~ 0. 

The process 7p — > J/0 X can take place through various mechanisms. The elastic 
process corresponds to X = p for which there are various models proposed, in particular the 
soft pomeron picture of Donnachie and Landshoff 0] and the Ryskin model |1| . The Ryskin 
model is based on perturbative QCD and is an attempt to describe a hard pomeron through 
a gluon ladder of the BFKL kind. The elastic process can be with proton dissociation or 
without. In both cases, with small momentum transfer, the J /if) particle retains the full 
photon energy (z » 1 with z = E^ J /E 1 in the proton rest frame). J/0 production with 
proton dissociation also leads to z values close to one even though its strictly speaking, an 
inelastic process. 

However the principle inelastic process is that of photon-gluon fusion which has been 



modelled by Berger and Jones |T7| who calculated this in perturbative QCD due to the hard 
scale given by the mass of the charm quark. In the colour singlet model for photoproduc- 
tion of J/0 which is based on the Berger- Jones picture, the formation of the J/0 state is 
accompanied by the emission of a hard gluon. A comparison of the colour singlet model to 
data gave a discrepancy in absolute magnitude which was attributed to missing higher or- 
der calculations (also called the "K-factor"). Subsequently, several improvements including 
a complete NLO calculation [T^| have been performed and compared successfully with the 
data. The possibility of colour octet contributions have also been studied - for an overview, 
see, for example, ||19| . All these various processes are shown diagramatically in Figure 12. 
There is not enough space to show the data for J/-0 production. The interested reader is 
referred to for data and experimental details. We will only briefly describe the data. 
The data show a rise with W 1P which in the HERA regime can be modelled by o 1P cx W^ p . 
A fit to the HERA data gives 5 = 0.64 ± 0.13 (including systematic and statistical errors). 
Using also the ZEUS data and lower energies, gives 5 = 0.90 ± 0.06. 
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The energy dependence from the Donnachie-Landshoff || soft pomeron model corre- 
sponds to 5 = 0.32. The Donnachie-Landshoff prediction falls well below the data by a 
factor of 3 or so. The prediction of the QCD-based model of Ryskin and including higher 



order corrections [16] gives much better agreement in general. The Ryskin model, based 
as it is on a 2 gluon ladder, depends quadratically on the gluon distribution taken at the 
scale (Q 2 + w3,)/4 « 2AGeV 2 . The distribution at this scale can be parametrised by 
xg(x) oc x~ x . For values of x ~ rn^/W^ p ~ 1CT 3 the gluon distribution from MRSA' |^TJ 
corresponds to A ~ 0.2 which gives good agreement with data between 10 and 150 GeV. 
The GRV parametrisation || corresponds to A ~ 0.3 — 0.4 giving a curve with a somewhat 
steeper energy dependence than the data. 

The inelastic data sample is selected by requiring 0.45 < z < 0.90 to reject elastic 
(diffractive) and resolved photon events. The inputs are the charm mass m c = lAGeV, 
Am S = 300MeV and the renormalisation scale is chosen to be the factorisation scale \f2m c . 
Once again the gluon distributions at low x can be parametrised by xg(x) oc x~ x with A 



ranging from (MRSD0') |22|] to 0.4 (MRSG) [^TJ. Approximate agreement is seen, however 
the data cover a wide range in absolute normalisation. Since the NLO corrections are not 
fully under control for p 2 — > and z — > 1, a restricted kinematic sample with z < 0.8 and 
p 2 > lGeV 2 has been studied by experimentalists. There is improved agreement for all the 



gluon distributions but the sensitivity is reduced. Again, details may be found in P0 |. An 
application of the asymptotic form of the gluon distribution as obtained from the double 
scaling limit to both the inelastic and elastic regions may be found in |2^ . 

To summarise, once again, the agreement with data is very good with standard gluon 
parametrisations, provided a hard pomeron (as that built up of gluon ladders as in the 
Ryskin model) is used. This picture explains the diffractive/elastic events far better than 
the soft pomeron version of Donnachie and Landshoff [[3] . In the elastic case the Berger and 
Jones picture using the colour singlet model of J/^ production seems to give a satisfactory 
description of the data once NLO corrections are included and suitable cuts on z and pf 
are incorporated. 
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7 Conclusion 



The structure function F2 has now been measured in the kinematic range 10~ 6 < x < 0.85 
and 0.1 < Q 2 < 30, OOOGeV 2 . The rise of F 2 with decreasing x is well explained by conven- 
tional DGLAP evolution but also by non- conventional dynamics like BFKL. A less inclusive 
measurement is required to distinguish these different evolution equations. The large ra- 
pidity gap events and diffractive scattering show that the (generic) pomeron is differrent 
in different processes but is harder than the Donnachie-Landshoff soft pomeron. Exclusive 
vector meson production also indicates that both for inelastic and elastic/diffractive events, 
perturbative QCD based descriptions explain the data fairly well. LEP and HERA have 
also given us the low x regime of the partonic structure of the photon, even though lack 
of space has not allowed me to discuss these aspects. The next frontier is that of polarised 
D1S at low x. The HERMES experiment is a step in that direction. 

The large body of theoretical work on low x has given us a good understanding of QCD 
dynamics at large parton densities. As a result, parton densities are now available over 
a large range in Q 2 and x. The results from HERA have given us new ways to connect 
partonic and hadronic degrees of freedom and contributed to a major improvement in our 
understanding of the short distance structure of hadrons. Low x QCD is one of the few areas 
of particle physics where theorists have had to work hard to keep up with experimentalists. 

Due to constraints of space, I have left out many other aspects of HERA Physics, like 
photon structure functions, charged current interactions, some aspects of jets in photopro- 
duction etc. The HERA Workshop on proton, photon and pomeron structure |24[] discusses 
in detail many of the issues touched upon here as well as some of the other experimental 
results that have not been discussed here. 
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HI Preliminary 1994 




Figure 11: x P F^ function of xjp in bins of j3 and Q 2 as denoted in the figure. The 
full line is the fit of an incoherent sum of the IP and / meson trajectories. The dashed line 
corresponds to the IP contribution, (from fT3|j ) 
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Figure 12: J/ip production mechanisms: a) elastic J/ip production via pomeron exchange 
b) diffractive proton dissociation c) elastic J/ip production in perturbative QCD (Ryskin 
model, see text) d) photon gluon fusion model for inelastic J/ip production in colour singlet 
model. 
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